% Written by: Orhan Torul
% Istanbul, Turkey, May 2017
function Gvec = sgu(ppp,bara1,bara2,bara3,bara4)

global a11g a12g a21g a22g a31g a32g a41g a42g a51g a52g a61g a62g a71g a72g a81g a82g;
global a11h a12h a21h a22h;

a11_g=ppp(1);
a12_g=ppp(2);
a21_g=ppp(3);
a22_g=ppp(4);
a31_g=ppp(5);
a32_g=ppp(6);
a41_g=ppp(7);
a42_g=ppp(8);
a51_g=ppp(9);
a52_g=ppp(10);
a61_g=ppp(11);
a62_g=ppp(12);
a71_g=ppp(13);
a72_g=ppp(14);
a81_g=ppp(15);
a82_g=ppp(16);
a11_h=ppp(17);
a12_h=ppp(18);
a21_h=ppp(19);
a22_h=ppp(20);
f_y_1=bara1;
f_y=bara2;
f_x_1=bara3;
f_x=bara4;
ga=[a11_g a12_g; a21_g a22_g; a31_g a32_g; a41_g a42_g; a51_g a52_g; a61_g a62_g; a71_g a72_g; a81_g a82_g];
ha=[a11_h a12_h; a21_h a22_h];

G =f_y_1*ga*ha + f_y*ga + f_x_1*ha + f_x;
Gvec=G(:);

% 
% 
% function deneme=BU(ppp,bara1,bara2,bara3,bara4);
% syms ag11 ag12 ag21 ag22;
% syms ah11 ah12 ah21 ah22;
% ag11=ppp(1,1);
% ag12=ppp(2,1);
% ag21=ppp(3,1);
% ag22=ppp(4,1);
% ah11=ppp(5,1);
% ah12=ppp(6,1);
% ah21=ppp(7,1);
% ah22=ppp(8,1);
% ga=[ag11 ag12; ag21 ag22];
% ha=[ah11 ah12; ah21 ah22];
% gav=ga(:).';
% hav=ha(:).';
% ghav=[gav  hav].';
% fa_y_tp1=bara1;
% fa_y_t=bara2;
% fa_x_tp1=bara3;
% fa_x_t=bara4;
% intera=fa_y_tp1*ga*ha+fa_y_t*ga+fa_x_tp1*ha+fa_x_t;
% deneme=intera(:);